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Introduction . 

For a first approximation the air flow around the airship 
hull is assumed to obey the laws of a perfect (i.e. free from 
viscosity) incompressible fluid. The flo-7 is further assumed 
to be free fror. vortices (or rotational motion of the fluid). 

These assumptions lead to very great simplifications of the 
formulae used but necessarily imply an imperfect picture of the 
actual conditions. The value of the results depends therefore 
upon the magnitude of the forces produced bv the disturbances in 
the flow caused by viscosity with the consequent production of 
vortices in the fluid. If these are small in comparison with the 
forces due to the assumed irrotational perfect fluid flow the 
results will give a good picture of the actual conditions of an 
airship in flight. 

* Dr. Max M. Munk's theory of the aerodynamic forces on an air- 
ship hull is presented in N.A.C.A. Technical Notes Nos. 104, 105 
and 106. This pajer was prepared by Dr. L. B # Tuckerman, a mem"?-:; 
of the special committee on the examination of the Naval Airship 
ZR-1, as a part of the Committee's revert and as an interpreta- 
tion and discussion of Dr. Hunk's papers. 
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General. 

The motion of a body through the fluid is accompanied with 
kinetic energy not only of its own motion but also of the motion 
of the fluid which it pushes aside. Since the fluid is assumed 
to be free from viscosity this kinetic energy of the fluid motion 
is not dissipated but accompanies the body in its motion, being 
transferred from portion to portion of the fluid as the body moves 
through it. The body, therefore, in any steady motion is accompa- 
nied by a steady configuration of fluid flow which changes only 
when the motion cf the body changes. If the velocity of the body 
is increased in any proportion the velocity of all portions of the 
fluid is increased proportionately (provided the velocities are 
small in comparison with the velocity of sound in the fluid; this 
is true here since the fluid is assumed to be incompressible) and 
the kinetic energy of the accompanying fluid motion remains pro- 
portional to the kinetic energy of the body itself. 

If, however, the character of the motion of the body changes, 
the shape of the accompanying fluid motion changes and the corre- 
sponding additional kinetic energy changes, although the velocity 
remain the same. 

Pure Translation. 

For a motion of pure translation Kirchhoff has shown that the 
kinetic energy (E f ) of the fluid can be written 

2 E f = p K x v£ + p KyVy + P K Z V| (1) 
where x, y, and z are three special axes in the body, mutually 
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perpendicular V x , V y and V 2 the corresponding components of the 
velocity of the configuration and p X x , p K y and p K z are 
"added inertias " corresponding to these three directions. On 
K x , K y and K 2 depend only the configuration of the body. The 
total kinetic energy E of the motion of the body is 

2E - 3E f + 2E b = ( pK x+ ni)Vx + ( OK y+T rJVy + (pK z+ m)V z 8 (3) 

Since no energy is dissipated, any change in the total kinetic 
energy of the motion of the body must be due to work done on the 
body (or by the body) 

- 6^ 6E=( pK x +Tii)V x 6V x+ ( pKyWVy 6V y + ( pX s +m)V 2 6V Z (3) 

If this change be due to a rotation, of the body without change of 
total velocity 

Vx + Vy + Vg - V" =' constant 
and V x 6V X + V y 6V y + V 2 6V Z - 0 

then - 6W = 6E = { p K x+ \)V X 6V x +( p K y+ \)V y 6V y+ ( pK z - \)V 2 6V Z U) 

-vhere the Lagrangean multiplier \ may be given any value -.re please. 
In order that there be no moment acting on the body tending to pro- 
duce this change it is necessary that 6 E - T 6 9=0 where 
T = the moment of force acting on the body and 5 0 the angle of 
rotation. This equation can obviously be satisfied (provided 
K" x j/b K ^ K 2 £ K x ) in three and only three mys. 
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K = - 0 K. 




0 



(5) 



X. = _ 



P K 2 , V x - V. 



V 



0 



These three mutually perpendicular directions in the body are there- 
fore directions of steady translation without the action of exter- 
nal moments. 

Lateral transfer of momentum . 

Consider a configuration of fluid flow A, (Fig. l) having 
a resultant momentum M in the y direction and no resultant mo- 
ment of momentum about the z-axis. Let this fluid mot ion be de- 
stroyed and replaced by an identical configuration in A r dis- 
placed a distance d having a component d sin 8 ("/here ? is the 
angle between the displacement and the direction of the i-os^ltont 
momentum) in the direction of the x-axis* To effect this cnari^e 
a negative resultant impulse -M must be applied to the fluid in 

and a positive resultant, impulse +M to the f luid in A 3 . 
That is, a resultant impulse moment Md sin 6 must act on the 
fluid. If, instead of a sudden transfer of momentum the transfer 
takes place continuously during the time t with a unif 
ity V such that d = Vt the impulse moment Md sin G is iue 
to a moment. 



acting during the time t. 

The distinction here between the momentum of the configuration 
of fluid flow and the momentum of a solid body should be noticed. 



T - 'IV sin e 



(6) 



- 5 - 



In a solid body the resultant momentum necessarily lies in the di- 
rection of its motion. The direction of resultant momentum of a 
configuration of fluid flow dce3 not necessarily coincide with the 
direction of the motion of the configuration. 

If T = 0 then 9 = 0 and the resultant momentum coincides 
in direction with the velocity. 

In the three mutually perpendicular directions considered 
above, since there is no resultant moment of force, the resultant 
momentum of the fluid must coincide in direction with the velocity. 
In these three directions therefore, the momentum of the fluid is 
given by 

M x - P K x V Xj M y = p K y V y , M z = P K s V z (?) 

and the resultant momentum in any other uniform translation is the 
resultant of these three moments. In general, the resultant mo- 
mentum M does not coincide in direction with the velocity of the 
body and thus needs a resultant moment T = MV sin G to be applied 
to the body in order to maintain a uniform motion of translation. 
This moment can be calculated either by 



JLE: ~a2.£ _ rn ( 6 ) 

(as in 4) or from T - MV sin G where M sin 6 is the transverse 
component of the momentum, (as in 6)- 

The calculation of the coefficients K x , K y and K z for any 
given body solves therefore for that body the problem of the total 
moments necessary to maintain it in uniform translation at any 
angle of pitch and yaw. 
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If the motion of the body is confined to the xy plane and 
K y = K 2 and K x K, , then 

3E f = DKj V* + 0K 3 V* = 0 (K s cos 3 a + K 2 sin* a)/ 
-here a is the angle of attack. Then 

T = -||f = 1/3 p V 2 sin 3 a (K 2 - K, ) (S ) 

or, other-rise, from equation (7) 

M x = p Kj. V x = pKi V cos a 
My = p K 2 V y - p K 2 V sin a 

and the lateral component of the momentum 

M sin a = !' v cos a - H x sin a 

- 1/2 p V sin 3 a (X 2 - X a ) 

and consequently, as before 

T = VM sin a = 1/3 p V 3 sin 2 a (K 8 - K a ) (9 ) 

Force Distribution . 

The determination of the force distribution which produces 
these moments requires a more detailed investigation. 
General Method . 

The general method may be sketched as follows: 
Under the assumptions here made the fluid flow possesses a ve- 
locity potential cp such that the component velocities of the fluid 
(not of the configuration) at any point are given by: 
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hiring determined this velocity potential the pressure at each 
point of the surface is evaluated from the extended Berrouii.li 
theorem 

p 1/2 v 2 -fi v = / v x J + v y 2 + v 2 J (10) 

Here Q is the. potential of the external forces acting on the flu- 
id. Since we are neglecting the change of pressure with height 
this may be treated as a constant. As Dr. Munk has shown, the 

£,Cp 

term 1ft may ' if desired * be transformed into 

— = _ v v cos e (11) 

where V is the velocity of the configuration at the point and 
9 the angle between the velocity of the configuration and;/the 
velocity of the fluid. 

This pressure is then integrated over the surface of success- 
ive zones of the ship, giving the resultant distribution of longi- 
tudinal and lateral forces along the ship. 

This process although perfectly general in theory is gener- 
ally impractical, since the velocity potential sP and consequently 
the velocity distribution has only been determined for a very few 
simple geometrical shapes, and even in these cases the computa- 
tions are laborious. 

Dr. Munk has, however, used the knowledge of the detailed 
pressure distribution based upon known velocity potentials in dis- 
cussing the effect of changing shape upon flow around two dimen- 



3 io rial shapes (No. 104, pp. 7,8 and 3\ 
A:-J ■ r o xiroate Soluti on. 

To avoid these difficulties, Dr. Munk attacks the problem ir. 
tnc following approximate way: The flow about any portion of the 
elongated ship is considered to approximate at any given instant 
the corresponding flow about an infinite cylinder having the same 
cress-section. (Fig. 3). In thia case the transverse added inertia 
is readily calculated from the well known case of two dimensional 
flow about an elliptic cylinder. 

The velocity potential in this case is determined from the 
complex function 

z = (x + i y) = f (w) = f ( <¥> + i t) 

where <p is the velocity potential and y the stream function. 
Here 

z = A w + - , 
w 

Proper choice of the constants A and B fits this to any elliptic 
cylinder between the limits of the infinitely thin flat plate and 
the circle. (See Lamb's Hydrodynamics, 4th edition, p. ?9, Lorenz.. 
Technische Hydro-Mechanik , p. 287.) 

If a and b are the major and minor semi-axes of the el- 
lipse, the added inertia per unit length p = p tt b 2 and 
p = p n a 2 . In the special case of a circular cylinder to 

which he confines himself in this presentation 

~ K<a = ^b = ^ J ~ S (12) 
-here S is the cross-section of the ship at this point, 



is of course zero. The contribution of any element of length dx 
to the total moment of the ship is therefore approximately from 

equation (9) 

£~ ±:=dT=l/2 P V s sin 2 c (Kg - K-J ) dx=l/ 3 p V 2 sin 2as dx (13) 

dT d 
since — = shear and -~ - lateral load per unit length, the 

^ x dx 
total moment T = 1/3 pV 2 sin 2 a/ S dx-1/2 pV 2 Q sin 2 o (14) 

where Q is the volume of the ship, and the lateral load 
F = / f dx is distributed according to the la-7 

f dx = ^4dx - 1/2 p V 2 sin 2 a ^ dx (15) 
dx dx 

This same method of reasoning he applies later to the problem of 
the rotating ship. 

The same result is arrived at more directly as Dr. Vunk ex- 
plained verbally, as follows: 

The transverse momentum of an element of length of the ship 
is, from equations (7) and (12) (Fig. 2) 

~ dx = dM - p V sin a S dx (16) 

dS 

If the cross-section S were increasing at the rate — the tran 
verse momentum rcould be increasing at the rate 

= PV sin a § dx = f dx. 
it at 

requiring a transverse load distribution f dx to impart this in- 
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crease of momentum. The equivalent of this increase of cross- 
section is imparted to the transverse air flow by the longitudinal 
component of the ship's motion (Fig. 3). As shown in the diagram 
the air which was flowing about the section S is after a time dt 

flowing about the section S + ^ dt where 4^ = V cos a 

dt at dx 

The corresponding increase of transverse momentum must be imparted 

to it by a laterally distributed force on the ship. 



f dx = p V sin a V cos a ~ dx 

dx 



- 1/2 p V 2 sin 2 a || dx 



(15) 



as before. 

The total moment on the ship calculated by this approximation 



was 

r2 



T - 1/2 P V 2 Q sin 2 a (14) 

obviously here the volume replaces the coefficient (K 3 - K T ) or 
equation (9). 

These coefficients K 2 and K a have been calculated for a num- 
ber of simple shapes. In particular, Lamb has calculated their 
value for ovary ellipsoids of different ratios of length to diameter. 

In this case, for all finite lengths K 2 - K 1 is less than 

the volume. Dr. Munk therefore proposes to apply a correction fac- 
te K 

tor (k 2 - k 3 ) (where k 2 = ^ and k 3 = ~) to the preceding for- 
mula, thus giving 

Total moment T = 1/2 p V 2 (k 2 - k a )sin 2 a Q (i?) 

Shear f~ - 1/2 P V s (k a - kjsin 2 a s (18) 
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Lateral force f dx = 1/2 p V 2 (k 2 - k,)sin 2 a g dx (19) 

where k_, and k x are Lamb's coefficients for the ellipsoid corre- 
sponding to the ship as calculated by the formula 

£ (ellipsoid) = ( J\ ^ ) (ship) (20) 

ROTATION % 
General 

If a body be in uniform translation parallel to one of its 
principal directions (V) , (Fig. 4), the added momentum of the 
fluid -vill have the same direction. About any axis A' perpendic- 
ular to this direction there Till be in general a resultant moment 
of momentum of the added momentum. There Till, however, be a line 
BB' parallel to the direction of the velocity such that the result- 
ant moment of momentum about any perpendicular axis (A) through it 
is zero. A similar line exists for translation in each of the other 
two "principal directions". These three lines do not in general 
intersect in a point. In bodies possessing certain types of aero- 
dynamic symmetry, however, they intersect in a point C, the aero- 
dynamic center of the body. If the body possesses geometrical sym- 
metry this aerodynamic center lies on the planes or axes of symme- 
try. This aerodynamic center exists in airship hulls and will be 
used as the center of reference for points in the body. The axis 
of x will be laid through it in the "longitudinal" principal axis 
of the body, this axis being an axis of central symmetry. 

The ship (Fig. 5) is supposed to be turning with a uniform 
angular velocity ~ about a fixed azis 0 where V is the linear 
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velocity of the aerodynamic center. The accompanying velocity con- 
figuration has a steady shape and steady speed and consequently a 
constant added energy but turns with the ship about the fixed cen- 
ter 0, The constancy of the energy requires that the resultant 
of all the forces acting on the ship pass through the center 0 
since otherwise the forces would have a moment about this a>:is and 
consequently add (or subtract) energy. These forces may be re- 
solved into a radial (centripetal) air force F p necessary to bal- 
ance the centrifugal forbe of the ship and of the accompanying 
fluid and a tangential (inertial drag) force F T either positive 
or negative, which is added to the frictional drag (neglected- here). 
The radial forces pass through 0, but the tangential force F t 
considered as applied at the aerodynamic center requires an accom- 
panying moment F T R to displace the line of action to 0. 

For the purpose of determining these forces the motion may be 
resolved into two parts, a parallel translation along the path and 
a rotation with angular velocity ^ about the aerodynamic center. 
If the center of mass of the ship coincides with its aerodynamic 
center this latter motion will involve no resultant forces nor re- 
sultant moments and consequently the resultant forces are calcula- 
ble from the parallel translation alone. 

The total tangential momentum M T (Fig. 6) of the ship in 

parallel motion is composed of two parts, ' r T due to the mass m 

i 

of the body 

M T = P V m (21) 
and M T due to the added tangential inertia 
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M T = pV (K 2 sin 2 a + K, cog- a ) (22) 

v;hile the total radial momentum M r is the added radial momentum 
alone and is 

M r = 1/2 p V (K a - K,) sin 2 o. (23) 
then (see Fig. 6) 

M x = M x sin P + M. r cos 9 M T = M T + M T 

My = M x cos 6 - Mr sin G 

From these the radial and tangential forces necessary to maintain 
the motion are 

F x = = ( k t cos G - M T sin 0)4f •= I (Mr cos P- M T - sine ) 

at at i\ 

F y - -gf =---(M r sin e + M T cos G )||= - £ (M r sin f + M T cos 9 ) 
If e = 0 F x = F T and F y = F r 

Then F T = | M r = 1/2 P V 2 ^ (K 2 - K a ) sin 2a (34) 
This represents a drag vfnsn a is positive. 

And F r I Mr = - p V 2 1 m - p Y ? 1 (K 2 sin 3 a + Kj cos 2 o ) (25) 

which is a centrifugal force. 

This computation is of course exactly the same as the usual 
calculation of centrifugal force in rigid dynamics, the only differ- 
ence being the existence of a transverse momentum, which gives rise 
to the "centrifugal" drag force. This is a generalized centrifugal 
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feros in the Lagrangee.ii sense. 

The drag F T is wholly due to air forces acting on the ship 

but of the centrifugal force F r that part due to the mass of the 

ship P V 2 -i m involves no air forces, the added centrifugal force 
R 

p V 2 ^ (K s sin 2 a + K, cos 2 a ) however, is transmitted to the 

ship by air forces acting on it. 

The drag Ft considered applied at the aerodynamic center is 
accompanied by the moment F T R = 1/2 P V* (X 3 - X,) sin 2 c ^hich 
is the same as the unstable moment in rectilinear motion (equation 
(9) ). The maintenance of the motion demands therefore (Fig. 7) 
a resultant force F and a moment T in addition to the aerody- 
namic forces here discussed. The fins alone supply the transverse 
component F 1 and the moment T = F'a. 
Distribution of these forces . 

Dr. Munk calculates the distribution of these air forces by 
the first method used in the case of rectilinear motion. Here, 
however, it is necessary to bear in mind that because of the curva- 
ture of the path the effective angle of attack of successive ele- 
ments of the ship's length are different. 

These angles of attack may be calculated as follows (See Fig. 8). 

x = _R =_R 

sin 6 siri(a f + | f ) cos a< 

a'=a-p=a- arc sin ( § cos a') 
Then sin 3a' = sin 2 a cos 2' (arc sin g cos a' ) 

- cos 3^ sin 2 (arc sin ~ cos a.') 



If a and # are both small this reduces to 

sin 2 c ' = sin 2a - 2 | ^ 
Then each element of length dx contributes an element of moment 

fM dx = 1/2 p V 3 (sin 2a - 2 | ) S dx (2?) 

The first term is due to the translation alone and the second term 
to the added rotation combined with the translation. Dr. Munk cal- 
culates these terms separately but the reasoning is equivalent to 
that here given. The total amount is 

T = 1/2 p V 2 sin 2 a / S dx - p V s I / S x dx (28) 

The first term is the unstable moment of the translational motion, 
and the second term is zero since / S x dx is the static moment 
of the volume about the aerodynamic center, which on the assump- 
tions here made coincides with the center of volume. As before, 
this calculation gives a resultant moment somewhat larger than 
acts on a ship of finite length so that he introduces again the 
correction factor (k 2 - k 3 ) in the first term. 

This factor gives the correct resultant moment. Since the 
remaining terms have no resultant, nor resultant moment, there is 
no obvious correction factor. Dr. Munk uses here k a * as a cor- 
rection factor instead of (k a - kj). 

The force dis + T' bvt i ^n is th^ 

♦Note: The difference is not great and it is all a matter of judg- 
ment but Dr. 'lunk's reason for using a different correction factor 
here is not clear to me. The forces are all calculated on the aam- 
basis of approximation. L. B. T. 
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~T d - f dx=l/3 P v'(k.-]f, )sin 2o gdx - P^(*ff + C)dx (29) 
and the total trans V3rce force 

F - / x dx - 0 

This approximate distribution of transverse air forces ther • 
fore accounts for the resultant unstable moment of the ship. 
It of course does not account for the drag. The undermined drag 
forces are, however, small, and being longitudinal, give rise to 
no appreciable bending moments in the hull. 

In addition, however, the approximation has yet to account fc 
the added centrifugal force (equation (25> ). 

P V 2 ~ ( Y s sin* a + ^ cos 2 a ) 

This force is of course snail since a and K 3 are. both small. 
For an ~ ratio of 6 and an angle of attack of 8 degrees it is 
less than 6 per cent of the ship's o^n centrifugal force. 

Of the two parts of this added centrifugal force, the first 

K v 
P V 2 sin 2 a = P V 2 — Q sin 2 a 
R R 

being due to the transverse added inertia can reasonably be assum- 
ed to be distributed according to the cross sectional area or 

f dx = p V 2 f§. S sin 2 a (30) 

The second term p V 2 — cos 2 a = p V 2 — Q 

R R 

(31) 



( a being small cos 2 a = 1 approx, ) 



requires a more detailed treatment, since its longitudinal distri- 
bution might give rise to considerable bending moments. As this 
term arises from the longitudinal added inertia alone, he consid- 
ers a case of longitudinal flow only, the flow arising from a sin- 
gle source and equal sink (Fig. 9). He chooses this flow (which 
gives a blunter airship model) instead of the corresponding ellip- 
soid because of the simpler mathematical treatment. The ccrre-^ 
sponding velocity is 

q> = f— - — ^1 (see Fig. 9) with the velocity distribution 

89 VD 2 / x - c x 4 c ^ 
x ~ ~ Ti 16~ V r a 3 " 4 t 

Here L = 3c + -i= approximately 

or nearly L = 2c. 

As may be seen from the indicated line of flow the longitudinal 
component of the velocity and consequently the added inertia is 
positive near the two ends but negative along nearly the whole of 
the side of the ship. At mid-section this negative velocity is 
approximately ^ diminishing to about —p- opposite the two 
sources and then rapidly changing sign around the nose. To simpli- 
fy his computation Dr. Munk assumes that it maintains its mid- 
section value ~ along the whole length and that the transverse 

2L 2 

velocity is negligible. This obviously results in an over-estima- 
tion of the bending moments produced. This flow, however, repre- 
sents a pure translation. The ship actually is rotating about a 



center 0 (Fig. 10), so that if V is the ship's velocity at the 
aerodynamic center, the surface velocity of the ship changes across 
the ship having a velocity V» = (V + y ^) at any point a horizon- 
tal distance y from the center. Dr. Munk* assumes that the air 
velocity retrains the same in the circular flight as in straight 
flight,** which gives 

* In a personal conversation, Dr. ISvnk states that this method of 
reasoning is different from the one he used, but as it arrives at 
the same result, is presumably equivalent to it. 

**?Jote: If the alternative assumption be made that the air velocity 
at any point of the surface in circular flight bears the came ratio 
to the surface velocity of the ship as it does in straight flight 
then 

. VD S ( . y> 

V * V (1 + 



) 



and p = _2J_!pi (i . -4-) (1 + if)* 

and the pressure gradient 



i*B. ,py 8 P! f 1 „-*>i d + 2>, ap r rcx^ 
dy R L* ' 4L^ / v R / R L ? 



imately. 



This pressure gradient is twice as great as on Dr, Munk's- assumption. 
It seems probable that the actual air velocity will lie between 
these two extremes > so that Dr. Hunk's assumption represents an un- 
der-estimation of the pressure gradient and consequently an under- 
estimation of the bending moments. As noted above, the assumption 
that the air velocity maintained its mid-section velocity 

— sr along the whole length, caused an over-estimation of the bending 
3L 

moments. These two factors rill cf course partially compensate each 

other, so that the bunk's assumption is probably more nearly correct. 



air velocity V - v - - -^pr 

configuration velocity V 1 = V (l + R ). 

Since the transverse air velocity is considered negligible 

9 - C and cos 9=1, then the pressure 'equation 
ooribined v;ith (11) ) 

p = - o v ' - oVI v> 008 6 + constant 



gives 



p = - - P - (iZr) + p [ 7 (i + ^] [ VDJ] + con3tant 

3 3L •' ' IV 2L 

and the pressure gradient". 

dy 3 It L 

This pressure gradient acts in the same -vay as a gravitatiu..: i 

0 V*"' D*^ 

pressure gradient due to a fluid of density — — t in a field of 

2 R L p^rf r 

horizontal intensity 1. The total lateral force is then ^ R L = 

(32) and is distributed along the ship proportional to the cross - 
sectional area. The added centrifugal force -t-— Qk, consists 
therefore of tvTo centrifugal forces §■ — Q (k, + gpr) ^ 33 ^ ccn ~ 
centrated practically at the ends of the ship combined with a can- 

tripetal force Q ~z (33) distributed along the ship projor- 

P v 2 

tional to the cross-sectional area. The factor — — Q is of cour 

R 

the centrifugal force of the ship itself, when in a state of static 

T D'" 
equilibrium. For k = 6, k 3 = .045 and -~ - .014. 

L> 3L 
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Transverse force on the fjns . 

A snail partof the centrifugal force can be balanced by the 
lateral viscous drag of the ship but the larger portion must be 
balanced by the lateral force en the fins. In addition, this later- 
al force must neutralize the unstable moment of the ship (Fig. 7). 
In his computation Dr. "lunk assumes this lateral force equal to the 
centrifugal force of the ship alone. This either neglects the add- 
ed centrifugal force or considers it neutralised by the lateral 
viscous drag. Equating moments (see Fig. 7) 

p V 2 Q | = p V 2 Q 1/2 (k a - k, ) sin 2 a 
or (kg _ kl ) sin 2a = j& (34) 

Summary . 

The lateral forces acting on the ship are then: 
1. The forces producing the unstable moment due to angle of 
at tack 



P V 2 



T = 1/2 p V 2 (k a - kj ) sin 2 a Q (17) 



a Q (35) 

The forces producing this moment are distributed according to the 
law 

f dx = a g dx (36) 

2. The lateral frrces due to rotation combined with tangential 
velocity. These forces have no resultant and no resultant moment. 
They are distributed according to the law 
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f dx * - fc, (x g + S) dx (39) 

5. The centrifugal forces or. the ship itself — — Q (35) 
provided the ship is in static equilibrium. If in addition the 
r.ass of the ship is distributed longitudinally- proportional 
to the cross section these are distributed according to the law 

,3 



P ^ - Q (31) 



f dx = -^y 1 - S dx \o<) 

These nearly neutralize the second term of (2)*. 

4. The added centrifugal force due to the added longitudinal 
inertia 

— iq 

TMs is distributed approximately as a concentrated load 

R 2 Ul + 2L 2 / 
at each end and a load distributed according to the la*r 

5. The added centrifugal force due to the added transverse 
inertia 

^tQ sin 2 o (25) 
This is distributed according to the law 



f dx k a sin*a S dx - S ; , a 3 cix ( 3C) 



♦Note: For any other mass distribution it ^rould be of course easy 
to calculate the corresponding force distribution. Since normally 
the static bending moments of the hull are everywhere hogging mo- 
ments, the actual force distribution is somewhat greater at the 
ends and less in the middle. L.B.T. 



6. The lateral force on the fins practically concentrated at 
the center of pressure of the fins 

- ^ Q (35) 

The sum total of all forces is then: 

Three concentrated loads 

a) at front end outward —r— ~t (k 1 + — ) (33) 

R 3 " 2 b 

oj at center of pressure — g — 0 (35) 
of fins inward 

c) at rear end outward —~ 2. ( k + > (3o) 

R 3 2L ' 

And a force distributed along the ship, with the resultant outward 
intensity 

f - X" C(a - k 2 x) || + (1 - k 2 -- - k 2 sin 2 a) s3 

(39) 

= £_X! [r a - k x) ^ + (1 - S a " k ^ ql 

r L a ks xj dx + 11 - *»- 3L 2 / b + r o^tp s] 

Ivote: 

The method of reasoning used in these papers introduces dis- 
crepancies between the computed forces and the actual forces due 
to two things: 

1) The viscosity of the air is assumed to be zero with the 
consequent elimination of all viscous drag. 

These discrepancies in the present state of the theory can 
probably only be estimated by comparison with experiment. 

3) The transverse flow about any element of the ship is as- 
sumed to be the same as that about the corresponding portion of an 
infinite cylinder. This assumption is most accurate when jy || 



is small. It will represent most closely the conditions amidships 

( 1 dS > The largest discrepancies will occur near the blunt 
v D dx ./ * 

nose of the ship ( f§ = and the next largest near the tail, 
1 dS 

where j ^ is finite out large. 

Since even small discrepancies in forces near the ends may 
result in relatively large discrepancies in the bending moments 
on the ship, it would seem to be very desirable to have some com- 
parison of the results of this approximate method with an accurate 
computation of the forces on a shape approximating that of the 
airship. 

The theory of the potential flow about an ovary ellipsoid is 
so complete that it is possible (although tedious) to compute the 
actual force distribution along such a shape both for straight 
flight and steady turning. 

It would seem that the comparison of the results of such a 
computation with the results of the approximate analysis given 
above would be of value in indicating the magnitude of the dis- 
crepancies involved. 

L. B. Tucker man. 
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Supplementary Note No, 1. ^odifi cation of D r. Muni's formulae. 

T!r. C. Eurgess lias called my attention to the practical 
disadvantage of an approximate load distribution which is not in 
equilibrium. By neglecting the added centrifugal forces in the 
calculation of the lateral force on the fins Dr. Munk leaves an 
unbalanced outward force of ^ — Q (k, + k 2 sin 2 a)/ This mak3 - 
no appreciable difference in the resulting moments on the ship 
but is inconvenient in practical computation, since it prevents 
the check obtained by computing both '*ays along the hull. 

This may be avoided by using the total centrifugal force in 
calculating the fin load, i # e. 



™-Qa(l + k,+ k 2 sin 2 o ) = p V P Q 1/2 (k s - kjsin 2 c 



or 



(k s - k 3 )sin 3a = d 4 kl+ i^sin- o) 

Since a is small the second approximation of its value 

77ill be sufficiently close for a numerical check. Then the total 

forces on the ship become: 



a) at bov- outward ^ l/2(ki 



4- 



3L C / 



b) at center of pressure P V f, 



of fins inward R. 



Q ( 1 + k,+ k 2 sin a a) 



c) at stern outward Q l/3(k, + ~g> ) 

and a force distributed alorg the ship with the resultant cut- 
■-ard intensity: 

f = -pT" [|a(l + k 1 + k r sin p o x _ k ? x \~ + ( 1-^ +k r sin" a -~s)S] 
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ou^plementary Note Ho. 2. ri3_cre par,cy Pgtws en Dr. f lunk's Theor y 

and N.P . L. Eat inmates . 

In computations for the ZR-l, Mr. Burgess has noted some dis- 
crepancy between Dr. bunk's theory and N.P.L. Estimates based on 
model tests. He pointed out that it is at least partially explain- 
ed by th^ neglect in Dr. Hunk's theory of the lateral resultant 
force on the hull arising from viscosity. The N.P.L. results show 



a 



Forc e on hull 
Total force 



400 _ . „ 



9500 



f - Force on f ins 6300_ „, r 
Total force" ~ 9S00~ ' " ,b 



Total force" " 9600" 
The lateral force on the hull is thus over 1/3 the total force 
and would make a considerable change in the results. 

It seems that the following method might give a somewhat bet- 
ter approximation. Assume forces as indicated in the diagram. 



Resultant 
force on 

f ins-. V 2 




hull=p._l Qh 
r R 



J. hen equating moments 



j TT~ Q a (f + P h) = C V ? Q l/2(k,-k 3 )sin 2 a, 



.?here Dr. Munk found 



then 



(k 8 - k, )sin 2c c =| 
sin 2 a i 



^~ _ f 

or q _ sin <d ~ G 



sin 2 a c 

sin 2 a 



= f + 3h 



Burgess gives a, - 7° 12' ; a 0 = 3 C 45» 

348? 

Substituting J . 73007 " ' 546 = 
values . 354 

The lateral forces on the hull have then apparently a result- 
ant applied about half way between the center of buoyancy and the 
center of pressure of the fin. 

It would seem then that a recoraputation by Dr. Munk's method 
based on an angle of yaw of 7° 12- with the addition of some rea- 
sonable distribution of lateral forces on the hull with a resultant 
at 0.51a might give a still closer approximation to the actual 
forces in a steady turn. 



Supplementary Note No. 3. Approxi mate .l9JM)£*±J& T J^£&.lS Coef- 
ficient;?. 

In comparing airships of different fineness ratio the vari- 
ation of Lamb's coefficients kj, k 2 and k, - k 3 may not always 
be negligible, although this variation need not be accurately es- 
timated. For such cases it may be rorth while noting the linear 
approximations given on the accompanying figures. These cover the 
whole range with a maximum error of 5% of the volume or the range 
4 < ^ < od with a maximum error of 27c. 

L . - . 

It is of course obvious that in the range 4 < g- < » earsboxic 

approximations would give still closer values. For instance, in 

„ /P>Y' 4 

this range the approximation k 2 - k a = 1 - 1.53 j has a ■ 
maximum error of less than 0. 34. In vie\7 of the roughness of the 
other approximations involved the -accuracy gained is probably net 
■worth the extra labor. 

L . S . la 
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Approximate values of Lamb's coefficients for prolate spheroid, 

ki, k ? and k s - ki 
Between ^ = 4 and oc. Maximum error is 0. 03 of volume 
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Approximate values of Lamb's coefficients k v , k r and k ? - kj fcx 
prolate spheroid. Error less than 0.05 part of volume over -hole 
range 
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